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1. INTRODUCTION 

Overhead lines are important components in the distribution of electrical energy. According to 
statistics, a large percentage of electricity costs, and most costumers’ outage, occur in the distribution 
network. In fact, about 80% of the customers’ outage are due to a failure in the distribution network [1]. 
Distribution networks in urban and rural areas often consist of overhead feeders. Poles, cross arms, 
insulators, jumpers, conductors and transformers are overhead distribution network equipment, thus forming 
a complex and repairable system. Conditions of all these component should be monitored, because the 
condition of each component will influence the overall condition of an overhead line [2]. Distribution 
companies strive to reduce their operation costs for achieving greater profits and maintaining their 
competitiveness [3]. In this regard, it is of great importance to reduce maintenance and repair costs, which 
constitute a significant part of the costs of operating the distribution network. Outages can be divided into 
four major categories. First, the failures that result from aging and failure of equipment and are related to the 
health of equipment and their service life. The second category is caused by wrong design and external 
factors and not related to the life and health of the equipment. The third category is due to the impact of trees 
and because the conductors are bare and also lack of timely trimming of the trees. Finally, the fourth category 
is caused by other factors such as human error and unknown factors. This information indicates that 
maintenance can only prevent one-third of the outages. According to experts, more than half of the unknown 
outage are most likely due to external factors. Therefore, to prevent the occurrence of nearly half of outages, 
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the design of the network must be changed and retrofitted against external factors. Therefore, failure rate 
modeling based on the age of the equipment does not seem to be realistic and is therefore a perfect model that 
takes into account the effects of environmental and external factors. As well as reflect the effect of 
performing corrective/preventive maintenance and reinvestment actions in failure rate improvement. 

Various methods have been used for evaluating and optimizing the risk, reliability and maintenance 
scheduling in distribution networks. As demonstrated in [4] utilized dynamic scheduling for long term 
maintenance management [5], [6], decision tree [7], [8], cost-benefit ratio [9]-[11], markov method [12], 
monte carlo [13]-[15]. Are the methods that risk and reliability of a network have been evaluated by them. 
The shortest path problem (SPP) method has been used in this study. SSP has been used for solving many 
different problems in the power system [16], [17]. Five actions are proposedfor managing risk in distribution 
networks, one option is to make no action during the network scheduling period, the next action is 
reinvestment for changing the design of the network and three actions are proposed for reinvestment. 


2. RISK MANAGEMENT ACTIONS 

The purpose of this study is to investigate the implementation of risk management in real 
distribution network. Risk management means reducing or keeping the risk at the optimal level using 
available and possible tools. But what tools do we have for risk management? Maintenance is only one of the 
risk management action, other tools such as replacement and reinvestment are available, and the choice of 
each of these action in each interval depending on the goals and constraints of each action. In describing the 
objectives, the expected level of risk can be mentioned. This issue that how much we want to reduce the risk 
depends on the goals and the prospects. Among the limitations are budget constraints and human resources 
that are influencive in achieving goals because we may not achieve the desired level of risk given these 
constraints. 

Our available and proposed actions for managing risk in distribution networks are illustrated in 
Figure 1. These actions include five categories. One option is to make no action during the network 
scheduling period. As a result, no budget are incurred, but the network situation worsens and the failure rate 
increases, imposing additional costs on the distribution companies, including the cost of the failures and 
energy not supplied. Another action is network maintenance. This means that we reduce the risk of network 
outages without drastically changing the network design by taking preventive and corrective actions. 
Maintenance can include tree trimming, minor maintenance and major maintenance. Therefore, the cost and 
effectiveness of each level of maintenance vary. The next action is reinvestment for changing the design of 
the network. Sometimes no level of maintenance can fix some of the outage, so it is essential that the design 
of the network be changed. For this purpose, based on the experiences of different countries and references, 
three actions are proposed for reinvestment as shown in the Figure 1. 


Risk 
| Management 


No Actions Reinvestment Maintenance 
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Minor Maintenance 


Figure 1. Risk management actions 


3. RISK MANAGEMENT MODEL 

The proposed risk management model is that we are faced with a choice of several actions in each 
period (e.g. a year). These actions include several maintenance actions and several investment actions that are 
selected according to the constraints and goals of the asset managers. The SPP is to find the optimum path 
from node | to end node n, which can be formulated as shown in the form of integer programming [18]: 
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where Z is zone number. In the (2) and (3), we have: indices i, j,k: index of node (1, 2,..., n); parameters n: 
number of nodes and c, ;: transmission cost of are (i,j); decision variables (binary) x; ;: the link on an arc 


Gj). 


3.1. Risk formulation 

The factors that are defined under the risk index include the risk of human injury and breaking 
safety and/or environment laws and the risk arising from a customer power outage [19]. There are also other 
common categories in the electricity distribution system such as financial risk, safety risk, environmental 
risk, supply quality risk, reputation risk, vulnerability risk and law risk [20]. Some references also categorize 
risk factors as limited to financial, safety, network performance and environmental risk [21]. 

In analyzing and evaluating different risk reduction actions it should be considered, that performing 
on not performing of each action what effect does it have on the risk factors. In this paper, only financial 
risks are considered. One of the economic indices of analysis for long-term scheduling is net present value 
(NPV) that considering discount rate, transfers all annual costs to the present year. Since the objective 
function is financial, therefore the analysis can be made through net present value function as shown in (4): 


n 
Cypv = » C x (1+ DR)” 
t=1 (4) 


where DR is discount rate, the economic or financial risk itself includes reinvestment, renovation costs, 
operation, maintenance, repairs, damage costs, costumer outage cost such as energy not supplied (ENS) costs 
and revenue loss cost. Note, because many scheduled maintenance programs perform in a hot line, we 
neglect the costs of the planned outage. 


NZ T ND 
Minimize CF = y y >: (cn) x y(t) + CagG x XF(t)) 
Z=1 t=1 k=1 


+ (CF(t) x A,(O)) (1 + DR~) 5 


Where y(t) is binary variable related to investment actions, X(t) is binary variable related to maintenance 


actions, oe T(t) is total cost related to the kth maintenance action of zone z in interval t and C/ (t) is total 
cost related to the outage occurred on zone z in interval t. 


3.1.1. Reinvestment cost 

Three actions are proposed for reinvestment ABC conductor, cover conductor and covering 
accessories of pole. To calculate reinvestment costs, we need to know the amount of capital needed to run 
each action per kilometer. Therefore, the equation of investing in the distribution network are as shown 
in (6). 


INV,* = inv* x L, (6) 
where k=1: covering accessories of pole; 2: cover conductor; 3: ABC conductor. 
3.1.2. Energy not served (ENS) cost 

ENS means the amount of energy demand measured in MWh which is not supplied in a given zone 
and in a given time period due to insufficient resources to meet demand. ENS costs include budget and labor 


constraints. It can calculate as shown in (7) and (8). 
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3.1.3. Revenue loss cost 
The economic or financial risk itself includes revenue loss cost. Total revenue loss cost is calculated 


by multiplying the ENS cost per unit by the quantity sold. It can calculate as shown in (9). 
CRevz = Ce X ENS,” (9) 


3.1.4. Maintenance cost 

The term maintenance cost refers to any cost incurred by an individual or business to keep their 
assets in good working condition. These costs may be spent for the general maintenance of items like price of 
one working hour necessary for the kth maintenance action or they may be used for repairs. It is calculate as 
shown in (10): 


cEMT Ct) = ( (CREE, x Labor*MT) + chMcome)) x 1, ) as 


where k is | for minor maintenance, 2 for major maintenance or 3 for tree trimming. 


NMx 


CkMcomp — > ckcomp 
J 
4 (1) 


As demonstrated in (10) means that each maintenance action includes more minor items. As an example for 
minor maintenance, the cost of non-standard span correction, replacement of fittings and correction of 
defective tie conductor is calculated separately or for major maintenance, replacement cost of failed 
insulators and poles, and cost of failed equipment are calculated separately. 


3.1.5. Failure costs 

Failure costs include price of one working hour necessary for repairing fault occurred, number of 
working hours necessary for repairing the fault, unplanned ENS Cost per megawatt-hour for zone z and 
revenue lost due to unplanned ENS resulted from outage occurred on zone z. It is calculate as shown in (12). 


CF(t) = ((Chunor * Labor®) + CRO™P) + (CENs, + CRE2) x (1+ q)** (12) 


3.1.6. Constraints 

Constraints refer to total cost related to the kth maintenance action of zone z in interval t, total cost 
related to the outage occurred on zone z in interval t, capital Investment associated with the kth investment 
action of zone z in interval t, number of working hours necessary for the kth maintenance action, number of 
working hours available in interval t for tree trimming and number of working hours available in interval t for 
minor and major maintenance or repair. All constraints can calcutale as shown in (13)-(17). 


NZ ND 
(eo x XE(t) + (CF) x a,(0))) < Budget(t) 

Z=1k=1 (13) 
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As demonstrated in (13) states that the total cost per year should not exceed the budget for that year. 
Equation (14) also emphasizes the limitation that the net present value of investment costs for all zones 
should be less than the amount of initial capital allocated to these zones. Also, in (15) notes that in long-term 
horizon (T), no more than one type of investment should be made for each zone. Equation (16) states that the 
number of man-hours for major and minor maintenance and emergency repairs per year should not exceed 
Labor table» the number available per year. Likewise, in (17) the rate of man-hours for tree trimming per 
year should not exceed the number Laborgiqitapie- 


4. STATISTICAL ANALYSIS OF OUTAGE DATA 

Statistical data in distribution networks provides operators, planners, and asset managers with 
valuable information on the performance and condition of the network and equipment. This data is collected 
in a variety of ways from a combination of software data such as computerized maintenance management 
system (CMMS), geographic information system (GIS), and communication networks such as supervisory 
control and data acquisition (SCADA) [22]. This information helps us to find equipment reliability indices, 
worst failure type, frequent failures, geolocation of failures and many more. 


5. CONDITION INDEX QUANTIFICATION 

Technically the best way to do this is to continuous monitoring of equipment condition. But in 
addition to high costs, this method requires a lot of infrastructure and effort [22]. In many cases the health 
index is defined to determine the condition of an equipment. The health index includes all the factors that are 
effective in determining equipment condition [21]. 

As illustrated in Table 1 the criteria are decoupled into four categories of vulnerability index, tree 
condition index, minor health index, and major health index. This decoupling has been made due to the static 
nature and dynamics of some indices over time, as well as to better reflecting the impact of each of the risk 
management actions on overhead line condition. 


Table 1. Check list for inspection of one zone of overhead lines 


Criteria Item Item Weight Item Score Criteria Weight 
Vulnerability 1 Vulnerability of poles to cars accident Wy S44 WwW, 
Index Vulnerability of conductors and pole accessories to birds Wy S12 

3. Vulnerability of conductors and pole accessories to other Wy3 S13 
external things 
4 Vulnerability of zone to lightning, Wy4 S14 
Weighted Mean of Scores bee Wij X S1j Xy =W, xX, 
X,= M 
dint Wij 
Tree 1 Vulnerability of conductors and pole accessories to trees Wo1 S24 WwW, 
Condition Weighted Mean of Scores YS Wo; X Spj peas 
Index x)= ye We =W,xX, 
jai 2) 
Minor 1 Condition of tie wires on insulators W31 S34 W; 
Health Index 2 Condition of connectors and jumpers W32 S32 
3 Condition of conductors distance W33 S33 
4 Condition of span sag W34 S34 
Weighted mean of scores bine Ws; X S3j Xutminor 
X; re way = W,xX; 
j= 
Major 1 Condition of poles (aging, corrosion) War Say W, 
Health Index 2 Condition of insulators Wa S42 
3 Condition of cut-outs, arresters and cable terminations W43 S43 
4 Condition of cross-arms Wa S44 
. M 
Weighted mean of scores . Dies Waj X Saj Xumajor 
Xa =W,xX, 


Overall Condition Index 


é:, Xy ots X tree+X u1_minor +X uma jor 


4 
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Another point in determining the factors affecting the overhead lines condition is the vulnerability of 
the overhead lines to external factors that are among the factors affecting the rate of outage or failure of the 
lines. The considering of this index was adopted after observing the extensive outages caused by these 
external factors in the statistical analysis. Table 1 was designed to determine the condition of a zone of 
feeders. The basis of zone division is the existence of maneuver points and the distance between two 
switching device. To quantification of the condition through this table, first each index is given a weight. A 
score will be assigned by the inspector for each index of each feeder zone. A score of 0 indicates the best and 
a score of | indicates the worst. The overall condition of each zone is calculated by the weighted average: 


Soe Xy + Xtree(t) + Xp1 major (t) + Xu minor (t) 


i ja W; (18) 


in the (18), X; is the overall condition of each zone as a function of time, Xy is the condition index of 
vulnerability factor, Xy; minor (t) is the minor health index and Xy; major(t) is the major health index of 
overhead line’s component. This will be more explained in next subsection. 


5.1. Health index as function of time 

In many cases the health index is defined to determine the condition of an equipment. The health 
index includes all the factors that are effective in determining equipment condition. Equipment such as poles 
and cross arms, fittings and jumpers, tie conductors, cutouts, and switches are subject to deterioration and 
aging. The following equation determines the health index of overhead line as a function of time [23], [24]. 


Xp (t) = Xun (e8ur(e—to) ) (19) 


In this equation, X;,;(¢) is an index of minor or major health at a particular time t in the future. Xj), is health 
index at time tgand f, is a constant coefficient for scaling the equation. 


5.2. Trees condition index as function of time 

Trees are fast growing and can reach a critical point in about one to two years. If the tree branch 
approaches the bare lines, it is very likely that a failure will occur. In fact, the vulnerability of bare lines to 
trees is very high. Here, we want to estimate the vulnerability of lines to trees in each zone using a time- 
based model. If the condition index is between 0 (best) and 1 (worst) we assume that the condition will fall 
from 0.05 to 1 within 2 years. With this information, the unknown parameters of the (20) are obtained. 


Xtree(t) = X trees (ePeree(t-to)) (20) 


In the (20), Xtree(€) is the condition index of the trees in each zone at time Xj;¢¢,, is the current condition of 
the trees at time tyand f;,e. is the unknown variable of equation. 


5.3. Scoring method of items 

After determining the weight of each item, there is a need to set a framework for scoring the items. 
Since Table 1 is designed for the inspector, it should be as easy as possible for the inspector to assign a score 
to minimize errors and preferences in scoring. Therefore, for each item in Table 1, we will provide a selective 
framework for scoring. For example, the method of scoring an item is described next subsection. 


5.3.1. Span condition 

To determine the item's score, the inspector must specify the number of spans with a non-standard 
sag and length. There are an average of 200 poles in each zone. So, based on knowledge of utility expert, we 
give the worst score (number 1) to a zone that has more than 10 non-standard spans. The scoring framework 
for this index is given in Table 2. 


5.4. Failure rate model 

The failure rate is a key parameter in risk analysis of distribution network and is known as an index 
and measure of reliability of the equipment and is also very closely related to the condition of the equipment. 
The worst condition of an equipment will result in the highest failure rate in the equipment [25]. In general, it 
can be said the best failure rate model is a model that, besides accurately estimating the failure rate, can cover 
all the factors involved in the outage of feeder, as well as reflecting the impact of maintenance, replacement, 


Int J Artif Intell, Vol. 11, No. 1, March 2022: 1-12 


Int J Artif Intell ISSN: 2252-8938 m) 7 


and all overhead lines preventive actions. In [25], a practical model is presented that, unlike common 
modeling which are as a function of time and age of the equipment, it models the failure rate of each 
equipment based on its equipment condition. The model presented is an exponential function of the failure 
rate as a function of condition. Equation (21) illustrates this function [25]. 


A(x) = Ae®*+C  failure/year (21) 


Table 2. Span inspection check list for each zone 
Score (0 to 1) Non-standard span (span/zone) _ Item 


Span Condition x= 1 0.05 
2<x<3 0.2 
3<x<4 0.4 
4<x<6 0.7 
6<x<8 0.8 
8<x< 10 0.9 
x > 10 1 


6. CASE STUDY 

The studied network consists of 10 feeders (19 zone) in Mashhad, that there are different types of 
feeders among them in terms of load type and geographical area. The set of customers who have an outage 
after each failure of the feeder and separation of the faulty part is defined as a zone. Therefore, a zone is 
considered between the two automation switch. 


6.1. Outage data analysis of case study 

Mashhad electricity distribution company (MCDC) with yearly distribution of 6661094 MWh, is 
one of the leading and important distribution companies in the electricity industry of Iran. However, there are 
various failures on the network each year due to different reasons. According to the information recorded in 
the information and outage management system of MCDC, about 4604 failures has occurred in the medium 
voltage feeders in Mashhad in the last five years, out of which 4190 were in the overhead lines. The all of 
these failures in the overhead medium voltage feeders Leading to 3816 MWh energy not supplied. 


6.2. Parameter calculation of failure rate model 

To calculate the parameters A, B, and C of (21) for each zone, we will follow the description of 
[25]. First, we select the initial values from the [25] based on the length of each zone. But since the outage 
rate of each zone depends on the length of the feeder and various factors, it is therefore necessary to calibrate 
these values by using the Chi-square error method [25]. 


6.3. Present condition calculation and condition model estimation 

To calculate the current condition of each equipment using Table 1, we obtain the weighted average 
of the characteristics and conditions of each zone in terms of its criteria and items. This average weight, 
which is the output of the table, indicates the condition of each zone under current conditions. The weight of 
the criteria and items is the same for its all feeders and related zones. To obtain minor and major £;,,from the 
(19), (22) and (23) is deduced respectively. Given the 20 year life span for the overhead lines and their initial 
state of 0.05, we obtain the By,. For to major and to,minor it is reasonable to assume that they are equal to the 
year in which repairs are done. 


Aig siinor 
© rae 
HIto0-minor 


Butmi = 
mee t — to minor (22) 
In (jmtar ) 
4 
Bur —_ HIt9—minor 
minor — 
ie t — tominor (23) 


Also to calculate Bi-e¢ by using the (20), (24) is concluded. We assume the worst-case condition Xtree = 1 
for the bare conductor as 2 years, for the covered conductor as 4 years and for the ABC cables as 5 years and 
the initial condition value for all feeders is 0.05. 
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x 
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Bree = (t a ie) (24) 
6.4. Effect analysis of risk reduction actions 

By performing each of the actions described above and applying the improvement effects of that 
actions on the condition index of the overhead lines in Table 2, the effectiveness of each actions is measured. 
There are now 5 actions outlined above, each applied separately and the rate of improvement in each 
condition index after applying actions is reported in Table 3. 


Table 3. Rate of reduction of decoupled condition indices after implementation of actions 


Risk Reduction Rate of Reduction of Decoupled Condition Indices After Implementation of Actions (%) 
Strategy Minor Health Index Major Health Index __ Trees Condition Index Vulnerability Index 
#1 (Action 1+2) 95 0 95 0 
#2 (Action 3) 0 95 0 0 
#3 (Action 1+2+4) 95 0 95 50 
#4 (Action 3+5) 95 95 95 70 
#5 (Action 3+6) 95 95 95 95 


Maintenance actions 

— Action |: trees trimming to prevent encountering branches with trees. 

— Action 2: minor repairs including replacement of damaged joints, correction of conductor spacing, 
correction of length and flash of span and tie conductors. 

— Action 3: general repairs including equipment repair, cross-arm replacement, defective insulators 
replacement or repairing damaged poles, and replacement aged equipment. 

— Reinvestment actions 

— Action 4: implementing covering of pole’s accessories: this action directly improves the network's 
vulnerability to birds, external objects and trees. 

— Action 5: replacing existing bare conductors with covered conductors: this action also has a direct 
impact on the amount of vulnerability to birds, external objects and trees. Of course, it is by itself 
vulnerable to storms and lightning. In the meantime and during implementation of the program, the 
connections, tie conductors, pole and cross-arms, the distance between the conductors and the bad span 
state will automatically corrected. 

— Action 6: replacement of current bare conductors with ABC cable: ABC cable not only eliminates 
conditions such as poor span condition, conductor distance, tie conductors, jumper and cross-arm, it is 
durable against external factors, especially trees, storms and lightning is and does not have the 
disadvantages of a covered conductor. 

Given that the overall condition of the feeder is resulted from the weighted average of the separated 
criteria, we therefore define that each action will only affect the self-relevant criteria. In accordance with 
Table 3, we define the strategies derived from the combination of the above actions. Rate of reduction of 
decoupled condition indices after implementation of actions are shown in Table 3 based on test and error 
using condition quantification framework in Table 1. 


6.5. Implementation of proposed algorithm on the case study 

In this section, the propsed model has been coded in the MATLAB [26] and solved by using 
priority-based genetic algorithm (priGA) 18] on the 10 feeders (19 zones) of the mashhad power distribution 
network described in the previous section, in 5-year horizons. The shortest path problem optimization model 
for 5-year horizon is shown in Figure 2. As it can be seen, we have 6 choices each year, depending on the 
specification of each zone and long-term scheduling, one of them will chose by the proposed algorithm. The 
chromosome dimensions of priGA after decoding for all zones will be NZ XT, i.e. 19*5 for the 5-year 
horizon. 


6.6. Problem assumption 
Typically, in any research, there are assumptions to bring the situation closer to the real operating 
conditions. In the following, we will explain the assumptions we have made to solve this problem. 
— Major maintenance will be carried out with the implementation of the ABC cable and coated conductor. 
— As a result of implementation of ABC and covered conductor action, minor health index and tree 
condition index are also improved. In fact, we do not charge a separate cost to improve these indices. 
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— Vulnerability to trees becomes less after running ABC cable and coated conductor. Therefore, the index 
of tree condition later returns to their bad condition. 

— Along with the implementation of the investment action for covering of pole’s accessories, the tree 
trimming action and the minor maintenance action is done too. 

— If we use ABC cable, there is no need for any major repairs until the end of the period. 

— Weare only allowed to use one investment action for each zone in each period. 

— Partial repairs are also done with the trimming. 


\ 


ye, 
as 


we = @——_-@—_@—_-@ —_® 


Figure 2. Shortest path problem for 5-year horizon 


6.7. Simulation result in 5 year horizon 

Simulation of the shortest path problem for finding the optimal strategy over the 5 years, is done 
based on model shown in Figure 2. In this problem, we have 6 options to choose each year. One option is to 
do no action on each zone, and the other 5 options include 5 different actions to improve the condition of 
each zone. By choosing an option for a year, for the following year only options can be selected that have a 
branch with the previous year option. In the following, we will examine the problems with considering 
constraints such as budget and labor. 

There are always some constraints for solving optimization problems that ignoring them, will made 
the solutions of the problems unrealistic. In this case, there are constraints such as budget and labor, which 
are determined by their optimal values in unconstrained case. The values of these constraints are given in 
Table 4. After solving the problem by considering these constraints, the paths presented in Table 5 are 
obtained. priGA convergence curve and optimal scheduling of action implementation are shown in Figure 3 
and Table 5 respectively. 

In this case, for all zones, the implementation of the investment has been result. In addition, 
investment on covered conductors (action 5) is proposed for zone 4. Investment on ABC cable is not 
recommended for any zone because of the high capital cost. In the cost analysis, we find that the total cost in 
this case has increased by 9% due to the constraints. The results of Table 4 shows that the annual cost and 
labor costs for maintenance and tree trimming did not exceed the defined limits. As is clear from the results 
of the unconstrained solution, many of the results are more than the set constraints. In the unconstrained 
solution, we found that due to disregard of constraints, there was little maintenance during the fourth, 
seventh, and tenth years, but instead, labor and cost density is in other years that may exceed the capabilities 
and budget and labor capital which shows the importance of considering constraints. 


Table 4. Result of optimization problem in 5-year horizon 


Year _ 1 2, 3 4 5 
Yearly Result Constraint 
Maintenance & Failure Cost (K$) 800 754.7 688.55 700.1 600.9 602 
Investment (K$) 800 649.4 165.3 8.4 0 0 
Man-hour Labor for Maintenance 9000 5572 8760 5074 3418 3983 
Man-hour Labor for Tree Trimming 300 227 91 222 188 209 
Cost Function (K$) 3,257.8 
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Table 5. Risk Management action schedules of case study comprising 19 zones in 5-year horizon 
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Figure 3. priGA convergence curve for 5-year horizon case 


7. CONCLUSION 

In this study, strategies have been proposed to improve risk, such as tree trimming, minor and major 
maintenance, covering pole accessories, replacement of bare conductors with covered conductors and ABCs. 
The shortest path problem model has been used to implement these strategies in the 5-year horizons. The 
results of applying the proposed method to the study network indicate that by considering of annual budget, 
available labor and allocated initial capital, zones with high loads, long repair times, and high failure rates 
require to replacement of bare conductor by covered conductor and other zones, they need to covering pole 
accessories. On average for all zones, minor maintenance and tree trimming are also done every other year, 
that confirming the validity of the modeling. The rate of major maintenance is also obtained on average every 
four years, which is consistent with the engineering logic of repairs. Overall, the results show the accuracy 
and precision of the proposed method. 
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